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ABSTRACT 



Context. Detecting regular dips in the light curve of a star is an easy way to detect the presence of an orbiting planet. CoRoT is a 

Franco-European mission launched at the end of 2006, and one of its main objectives is to detect planetary systems using the transit 

method. 

Aims. In this paper, we present a new method for transit detection and determine the smallest detected planetary radius, assuming a 

parent star like the Sun. 

Methods. We simulated light curves with Poisson noise and stellar variability, for which data from the VIRG0/PM06 instrument on 

board SoHO were used. Transits were simulated using the Universal Transit Modeller software. Light curves were denoised by the 

mean of a low-pass and a high-pass filter The detection of periodic transits works on light curves folded at several trial periods with 

the particularity that no rebinning is performed after the folding. The best fit was obtain when all transits are overlayed, i.e when the 

data are folded at the right period. 

Results. Assuming a single data set lasting 150 d, transits from a planet with a radius down to 2 Re can be detected. The efficiency 

depends neither on the transit duration nor on the number of transits observed. Furthermore we simulated transits with periods close 

to 150 d in data sets containing three observations of 150d, separated by regular gaps with the same length. Again, planets with a 

radius down to 2 Re can be detected. 

Conclusions. Within the given range of parameters, the detection efficiency depends slightly on the apparent magnitude of the star 

but neither on the transit duration nor the number of transits. Furthermore, multiple observations might represent a solution for the 

CoRoT mission for detecting small planets when the orbital period is much longer than the duration of a single observation. 

Key words. Stars: planetary systems - Occultations ~ Techniques: photometric - Methods: data analysis 



1. Introduction 



Beginning with the discovery by iMayor & Oueloz I d 19951) 1 2 
years ago of a giant planet, 5 1 Peg b, around its Sun-like host 
star, the search for planets outside the solar system has become 
a major topic in astrophysics. Among the various methods avail- 
able to search for exoplanets, the transit method presents the op- 
portunity of directly determining the planet's radius relative to 
that of its parent star and, if more than one transit is observed, 
its orbital period. It also has the advantage of beeing possible to 
monitor many thousands of target stars simultaneously. This is 
of course necessary, since the probability that a planet transits 
its star depends on the radii of star and planet, as well as the or- 
bital distance and the eccentricity jRauer & Eriksson1l2007h . It 
is never a large number, viz. 



poorly known factor rjpianet (number of planets per star system), 
one realizes this need for (simultaneous) observations of large 
numbers of stars. 

Diminishing of the light curve relative to the light from the 
star, the photometric signal P, is of course immediately translat- 
able into the ratio of planetary to stellar radius 
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where a is the semi major axis of the planetary orbit, rj and rp 
are the stellar and planetary radii respectively, e the eccentricity, 
and CO the longitude of the periastron. It is easy to see that, for 
an Earth analogue, the probability of a transit beeing observed is 
about 0.5%, while for a similar planet at 0.05 AU, the geomet- 
rical probability rises to about 10%. When also folding in the 
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Therefore, the minimum size of the planet that can be detected 
will only depend on the precision with which one can measure 
the photometric output from the star and know the stellar radius. 
Successful results have been obtained for a score of tran- 
siting objects, since the first detection by ICharbonneau et all 
(l2000h . but all of these are essentially Jupiter-like objects in 
very close orbits (i.e. less than a few days). If we want to de- 
tect and study planets more akin to our own Earth, we are fac- 
ing formidable problems, where the stellar radius currently can 
only be discerned from model calculations, and even then is still 
contains severe uncertainties. The achievable photometric pre- 
cision is governed by the problem of acquiring long and un- 
interupted sequences of data. When observing from the ground 
we are faced with enormous challenges. In contrast, by launch- 
ing a telescope into space, into an orbit allowing these long, 
uninterrupted data acquisitions, much higher precision can be 
achieved. Nevertheless, by improving the photometric signal, 
one also increases the ambition. Thus, with the advent of space- 
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based assets, the goal is now set at the actual detection of a planet 
with the same radii as our own Earth. 

Simultaneously, during the development of the first space 
missions, a number of transit-detection algorithms have been 
proposed in the l i teratur e including Bayesian al gorithms (i.e. 



Ai grain & Favata I |2002|) . mat ched filters (i.e, Jenkins et alj 



19961) . box-shape transit finder ( Aigrain & Irwin '2Q04), and the 



box-fitting least squares method ( Kovacsetal. 2002;) . A com- 
par ison of t h e capab ilities of existing methods was first made 
by Tinglev (l2003al) (revised in Tingley 2003b) and later by 
iMoutou et al.l (l2005h on updated versions of the algorithms. 
Besides the detection method, filtering of pho ton noise and 
stellar variability also play an important role. In ICarpano et al.l 
(l2003h . we used an optimal filter to significantly reduce the stel- 
lar variability and enhance the S/N ratio of the transits. However, 
this filter uses a reference signal and hence some a priori infor- 
mation of the shape of the signal we are looking for. 

The CoRoT spacecraft was launched on December 27, 2006 
and has been carrying out routine scientific operations since 
February 2, 2007. CoRoT stands for COnvection, Rotation ,and 
planetary Transits. Briefly, this satellite is carrying out ultra 
high-precision photometry for the dual purpose of studying as- 
teroseismological signals from stars and both discovering and 
studying extra-solar planets with unprecedented precision. The 
focal plane is divided into two segments, one each for these ob- 
jectives. Each of the parts of the focal planet (designated 'exo' 
and 'seismo') is equipped with two CCD-based detectors. The 
exoplanetary detector records the light curve of up to 12 000 tar- 
gets. Of these, 20% have a time sampling of 32 sec, while the 
rest are sampled every 512 sec. Continuous monitoring of up to 
150 days for each selected field is carried out, and a grism in 
front of the detector allows thousands of the light curves to be 
obtained in three colours for stars with visual magnitudes from 
about 12 to 16. For further basic information about the mission 
we refer to ESA SP-1306 (2006, The CoRoT Mission) for a de- 
scription of all aspects as of pre-launch. 

One of the main goals of CoRoT is to detect smaller exo- 
planets than previously possible with any method. Radial veloc- 
ity measurements is currently limited to the detection of plan- 
ets 5-10 times as m assive as the Ear th itself in orbits around 
very low-mass stars JUdrvet alll2007h . The occultation method 
as carried out from the ground also has limitations imposed by 
both the atmosphere and, more important, by the interruptions 
inherent in any lightcurve obtained from the surface, although 
extensive networks developed in recent years alleviate this. In 
this case, the depth in the light curve that can be obtained from 
the ground u nder the best of circumsta nces is at best a few tenths 
of a percent JRauer & Erikssonll2007l) . representative of planets 
of either the Jupiter or Saturn class. In the case of CoRoT, on 
the other hand, it has been shown by in-orbit verification that the 
data is photon-noise-limited, or almost so, in the whole magni- 
tude range for the exo-planetary segm ent for integration times 
shorter than 1 hour dBaglin et alll2007l) . The detections of Earth- 
sized planets then becomes possible with the CoRoT mission - 
at least for shorter orbital periods - but is nevertheless is a chal- 
lenging task. 

In this paper we present a new, although simple, transit de- 
tection method that could be applied to CoRoT data after hav- 
ing applied a smoothing and detreding filter. We do not attempt 
to compare the efliciency of the method with respect to other 
ones, for which a proper analysis should be performed much like 
the one in Tinglev (2003a) and Moutou et al. (2005J). Here we 
concentrate on describing the method, and quote its capabilities 
and its efficiency at detecting small planets when few transits 



are observed. The detection is applied to simulated light curves 
with Sun-like stellar variability, lasting 150 days as is the case 
for CoRoT. Assuming the radius of the parent star is fixed to 
solar value, we check the radius the planet needs for beeing de- 
tected. Stellar apparent magnitude, transit duration and orbital 
period, hence the number of transits, are all taken to be random 
in a given range of values. To avoid a long computation time, we 
limited our simulation to 100 light curves for each stel lar activ- 
ity leve l and therefore did not follow the approach of lTinglevI 
(l2003al) . which uses Monte Carlo simulations. 

We also tested the transit detection method in the case of 
three observations of the same object separated by 150 days, 
where only one transit is visible in each data set. This simulates 
the observation of the same region at several separated epochs. 

The paper is organised as follows. Section [2] describes the 
way the simulated input light curves were derived. In Section [3] 
we describe the filtering and the transit detection method we 
have applied. The performances of the method are evaluated 
in Sects. |4] and |5] for the single and multiple observations, re- 
spectively. Finally, the discussion of our results is presented in 
Section |6] 

2. Data simulations 

To test the efficiency of our transit detection method in finding 
Earth-sized planets around a Sun-like star, we simul ated light 
curves using VIRGO/PM06 data (Frohlich et a l1ll997h from the 
SOHO telescope as done in ICarpano et al.l (|2003|) . Since data 
from VIRGO are only available on a hourly time grid, we inter- 
polated them to have a data sampling of lOmin. The total du- 
ration of the light curve is 150 d as expected from the CoRoT 
mission (where the time resolution is 32 s or 512 s). To simu- 
late a high level of stellar activity, VIRGO data from the year 
2000 were used, while low activity was represented by data from 
the year 2006. The corresponding VIRG0/PM06 data sets are 
shown in Figure [1] On the top of these data, we added Poisson 
noise, assuming that AF/F -7x 10""^, where F is the total num- 
ber of photons collected in 1 h for a star of apparent magnitude 
15.5. This is what is expected for the CoRoT mission. No instru- 
mental noise of any kind has been introduced into the simulated 
light curves. 

Dips associated with the transit of a planet are simulated us- 
ing theUniversal Transit Modeller (UTM) software, developed 
by iDeegl (Il999l) . The utilised parameters values are listed in 
TablelU the mass, M*, the radius, R*, and limb-darkening coef- 
ficient, Ld, of the parent star, which are kept fixed to solar val- 
ues. The radius of the planet, Rp], the orbital period, P, and the 
phase, ph, (thus the number of transits), and the transit duration 
change randomly within a given range of values. We require that 
at least 3 transits are present in the light curves, with a maxi- 
mum of 6 transits. The apparent magnitude, ranging from 12 to 
16 is related to the Poisson noise added to the light curve. An 
example of the output of the UTM code, zoomed around the first 
transit, is shown in Figure|2]where Rp\-3 R®, P=37.5 d, Ph=Q.5, 
tdur-lOh. The sharp drop and increase in flux are caused by the 
ingress/egress of the planet, while the smoother variation close 
to the bottom is due to limb-darkening effects. 



3. Filtering and periodic transit detection 

3.1. Filtering 

To filter the Poisson noise and detrend the light curve for stel- 
lar variability, we used a low-pass and a high-pass filter, respec- 
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Fig. 1. VIRGO/PM06 data from the SoHO satellite, from year Fig. 3. Simulated light curve before and after filtering. Transits 
2000 (top) and 2006 (bottom) used to simulate stellar activity at are built with the same parameters as in Figure|2] stellar variabil- 
high and low level, respectively. ity is taken at the high level, and y= 12 mag. 



Table 1. Parameters used in the Universal Transit Modeller soft- 
ware to simulate the transit of a planet in front of its parent star. 



Parameter 



Value/Range of Value 



Star mass, M^ 1 Mq (fixed) 

Star radius, R^, 1 Rq (fixed) 
Limb darkening coefficient, Ld 0.51 (fixed) 

Apparent magnitude, V [12-16] mag 

Transit duration, tdur [4.5-36] h 

Planet radius, Rpi [1-5] Re 

Orbital period, P [25-50] d 

Orbital phase, P/i [0-1] 



1.0000 




Z 0.9994 



0.9992 - 



17.0 17.5 18.0 18.5 19.0 19.5 20.0 20.5 
Time (days) 

Fig. 2. Output of the UTIVl code, zoomed around the first transit 
, where 7?pi=3 Re, P=37.5d, Ph^O.5, tdur^lOh. 



lively, equivalent to a standard box-shaped filter in the Fourier 
domain. The disadvantage of these filters is that the shape of the 
dip is slightly changing. Without passing through the Fourier do- 
main, the data sets are filtered by a boxcar smooth over a certain 
width, using the IDL tool smooth. To reduce the Poisson noise, 
we smooth the data over a width of 25 (in units of 10 min). This 
width, linked to our data sampling, limits the lower value of our 



simulated transit duration (4.5 h). Outputs from the CoRoT mis- 
sion have a data sampling of 32 s (in the oversampled mode), 
therefore decreasing the lower limit for the transit duration to be 
detected. The smoothing unfortunately reduces the sharpness of 
the transit ingress/egress. This shape's characteristic can be used 
to distinguish dips from an eclipse with respect to dips from stel- 
lar variability; however, this will have no effect on our detection 
method as described below. 

The detrending is performed in two steps: we first smooth the 
data over a large width (wider than the maximal transit duration, 
i.e. 1.5xfrfMrn,ax=54h) and then remove the long-term smoothed 
light curve to the previous short-term smoothed one. This results 
in a light curve filtered for stellar variability and Poisson noise. 
Extreme values are set to the mean light curve value (1 if nor- 
malised). An example of a light curve before and after filtering 
is given in Figure|3] where the parameters for the UTM code are 
taken as in Figure |2] Stellar variability is assumed to be high, 
while the apparent magnitude V - 12 mag. 

3.2. Periodic transit detection metliod 

Our transit detection method has the advantage of beeing very 
simple. We first fold the light curve over a given range of trial 
periods. The only diff erence with respe ct to the conventional 
epoch-folding method dLeahy et al.lll983h is that the data are not 
rebinned after folding. On the folded data set, for each trial pe- 
riod, we first sort the data in an ascending order for the phase 
values and then fit a Gaussian function using the IDL function 
gaussfit. Results are shown in FigureH) when data are folded 
at the correct period, the Gaussian fits the data relatively well 
(top of Figure nil. In contrast when data are not folded at the cor- 
rect period, the width of the Gaussian function is constrained by 
the data sampling, which results in a bad fit. This is also the case 
when folding at half of the period, one third of it, etc. Testing 
over a full range of trial periods, the best value is given by the 
minima of the chi-square function as shown in the next section, 
in Figure |5] Note that the fit is performed over the entire folded 
light curve, while in Figure |4] only the region around the transit 
dip is shown. The shape of the transit is modified by the filter, 
which creates bumps before and after it. 
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Fig. 4. Light curve folded at the correct (top) and wrong (bottom) 
period, without any rebinning. Overlayed is the best fit Gaussian 
function (see text for more details). 



Other minima can be observed at multiple values of the cor- 
rect period. If the S/N ratio is not too low or if the period step 
value is not too high, the chi-square value associated with these 
sub-harmonics will always be at the secondary minima. No sig- 
nal will be found at half of the period, one third of it, etc., and 
thus no harmonics are present in the periodogram. The period 
step should be chosen as small as possible, ideally one time unit, 
or at least l/«th of the transit duration when no more than « tran- 
sits are expected. 

This detection method presents the disadvantage that the po- 
sition of the fitted Gaussian function is given by the minima of 
the filtered light curve. If this minima is not associated with any 
of the transits but with some dip coming from the stellar activity, 
for example, the correct period will not be found. 

Another property of the method is that the detection effi- 
ciency does not depend on the number of transits observed. The 
right period of a data set containing a large number of tran- 
sits will be detected in the same way that with only a few ob- 
served occultations. However, the presence of many transits will 
increase the confidence level of the period detection (see next 
section). 

Comparing with other methods qualitatively, this one has the 
advantage of not needing to create a set of test light curves to 
be compared to the observed one as do ne for the matched-filter 
jjenkinset al.lll996l) or the method of iProtopapas et alj ( l2005b 
(which uses analytical function to characterise a transit). It does 
not need any assumption on any parameter, like the transit dura- 
tion, as with t he sliding tran s it tem plate correlation developed 
by Team 1 in iMoutou et al.l (l2005h . The non-improvement of 
period detection with an increase in the number of transits is 
also o bserved in the method used by Team 2 in iMoutou et alj 
(l2005h . They make a box search for all data points deviating 
from the average signal by 3 sigma, remove spurious events, 
and subsequently search for a periodicity in all detected epochs. 
Our method has the advantage of doing everything at once. 
Furthermore, we do not detect harmonics of the periodic sig- 
nal but only the sub-harmonics, which decreases the noise level 
in the chi-square function. Of course a proper comparison would 
be necessary to test the efficiency of the periodic transit detection 
and/or time consumption of this method with respect to existing 
ones. 
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Fig. 5. Chi-square ratio values for the Gaussian fit of light curves 
folded over trial periods ranging from 5 h-55 d, for an input pe- 
riod at 50 d (top) and 25 d (bottom). The minimum of the curve 
gives the value for the best-fitted period. 



3.3. Confidence levels 

To determine the confidence level of a period detection, we esti- 
mated the ratio between the chi-square resulting from a linear fit 
and the one resulting from a Gaussian fit. We call this the y^ ,. 

^ to ^ ratio 

(see Figure|5] top, for three transits, and bottom, for six transits). 
The second chi-square increases with the number of transits, but 
the first increases faster Thus when increasing the number of 
transits, the v^ ,. function goes to lower minima, therefore en- 

' '^ ratio to ' 

hancing the confidence level. We then compared this minimum 
to that expected for transit free data. In practice, we simulated 
100 light curves without transits using solar data both at high 
and low activity levels. We filter the data, search for a periodic 
signal transit, and calculate the ;if^3(jg. It turns out that 68% of 
the light curves will provide ^X^-^iio -' 0.971, while 90% will be 
> 0.989 and 99% greater than 0.997. Thus a transit is detected at 
a confidence level of 90% if the ;if ^^^^^ > 0.989. These thresholds 
do not depend on the activity level of the star since they have 
been pre-filtered. 



3.4. Estimate of the planetary radius 

Since we kept the radius of the parent star equal to 1 Rq, the 
depth of the transit will directly estimate of the planet radius. 
This depth is given by the amplitude of the fitted Gaussian func- 
tion. In Figure |6] we show the correlation between the simulated 
planet's radius and the one determined by the Gaussian fit, when 
the right period has been recovered. The maximal deviation be- 
tween the simulated and detected radius is ~ IR®. 

We also tried to get an estimate of the transit duration, from 
the width of the fitted Gaussian function; but due to the high 
level of noise, no correlation was found for these two parameters. 
However, in a second step once the transits have been detected, 
a deeper analysis of the source allows a more accurate fit. An 
estimation of the transit duration is then provided by the width 
(2x4 sigma) of the Gaussian function. 
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Table 2. Performance of the periodic transit detection algorithm 
in case of high stellar activity. 

Planet radius (Rq) 



Radius Simulated 

Fig. 6. Correlation between the simulated planet's radius and the 
one determined by the Gaussian fit, when the right period has 
been recovered. 

3.5. How the method deals with unevenly sampled data or 
with short/long data gaps 

The detection method described above has the advantage of not 
being too affected by unevenly distributed data. The prior fil- 
tering of the light curve, on the other hand, works properly 
only for evenly sampled data. This problem can be solved eas- 
ily by rebinning the uneven data set on to a regularly-spaced 
grid. Concerning the gaps, if short ones are present (short with 
respect to the transit duration), this will again affect the results 
only in the pre-filter part. A solution will then consist in interpo- 
lating the data over the missing values. If medium-sized gaps are 
present, the procedure and results will stay the same as long as 
no transits are present in the gaps. If transits are present in these 
gaps, the correct period will not be found even if the dips are 
very pronounced. When long-lasting gaps are present, each data 
set should then be pre-filtered separately. The detection method 
works also well with long gaps, depending on their length and 
the orbital period to be detected. This is of course important if 
we want to search for planets with periods longer than x 1/3 
of the length of the uninterrupted data, and thus have to return 
to observing the same object after an interruption that could be 
up to the same length as the original data set (e.g. because of 
reasons of orbital mechanics with respect to the spacecraft). An 
application of this case is given in Section|5] 

4. Performances of the detection method 

To test the efficiency of the detection method, we simulated 100 
light curves as described in Section|2](Table[T]i. Results are dis- 
played in Figure [T] for the case of high stellar activity and in 
Figure [8] when stellar activity is low. We compare the minimum 
size of the detected planet as a function of the apparent magni- 
tude (top), the transit duration (middle), and the orbital period 
(bottom). Table |2] shows the total number of simulations (NB., 
1st column) for a given range of the input parameters of Table [1] 
radius of the planet, magnitude, transit duration, and period (re- 
lated to the number of transits), for high stellar activity level. The 
2nd, 3rd, 4th, and 5th columns show, respectively, the number 
(in %) of correct period detection (YES), of failed period detec- 
tion (NO), of false alarm (FA.), and of missed detection (M.D.). 
Table[3]shows the same results but for a low activity level. 



Radius 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


1-2 


22 


23 


77 


9 


13 


2-3 


27 


48 


52 


22 





3-4 


31 


94 


6 
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10 


4-5 


20 


100 











Magnitude (Mag) 


Mag. 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


12-13 


26 


81 


19 





4 


13-14 


24 


79 


21 


4 


4 


14-15 


25 


72 


28 


8 


8 


15-16 


25 


36 


64 


24 
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Transit duration (hours) 


Tdur 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


4-12 


22 


77 


23 


9 





12-20 


28 


68 


32 


11 


7 


20-28 


26 


65 


35 


8 


8 


28-36 


24 


58 


42 


8 
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Period (days) 


P 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


24-30.5 


28 


75 


25 


11 


7 


30.5-37 


26 


62 


38 


12 


4 


37-43.5 


25 


68 


32 


8 


4 


43.5-50 


21 


62 


38 


5 


10 



Table 3. Same as Table|2]for low stellar activity level. 
Planet radius (Rq) 



Radius 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


1-2 


24 


8 


92 


8 


4 


2-3 


25 


60 


40 


12 


12 


3-4 


30 


93 


7 








4-5 


21 


100 
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Magnitude (Mag) 


Mag. 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


12-13 


32 


75 


25 








13-14 


26 


65 


35 


4 


4 


14-15 


20 


80 


20 





15 


15-16 


22 


41 


59 


18 


5 


Transit duration (hours) 


Tdur 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


4-12 


21 


71 


29 


5 


10 


12-20 


27 


78 


22 





7 


20-28 


25 


56 


44 


8 





28-36 


27 


59 


41 


7 


4 


Period (days) 


P 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


24-30.5 


27 


59 


41 


7 


7 


30.5-37 


30 


60 


40 


7 


7 


37-43.5 


24 


67 


33 


4 





43.5-50 


19 


84 


16 





5 



The main results are that, when assuming solar variability 
and Poisson noise level as given in Section |2] planets with a 
radius down to ~2Re can be recovered when the parent star 
has a radius of 1 R©. The performance depends slightly on the 
apparent magnitude but not on the transit duration. Finally, as 
expected, the efficiency does not depend on the orbital period, 
which is directly correlated to the number of transits. 
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Fig. 7. Results of searching for small planets in 100 lightcurves. 
Performance of the periodic transit detection algorithm in case 
of high stellar activity. Dots are associated with a successful re- 
covery of the period with a confidence level > 90% and circles 
with a confidence level < 90%. A failed recovery of the sim- 
ulated period is shown by a cross when the confidence level is 
< 90% and by a star for a level > 90%. Parameters and units are 
explained in Section |2l 



5. Multiple observations with long and regular gaps 

In this section, we investigate one possible application for the 
CoRoT mission: what if a planet has an orbital period close to 
1 yr, when a transit is observed in the first 150 d data set. Since 
the satellite cannot observe the same field again before another 
150 d, would it still be possible to recover the orbital period of 
the planet? Since the detection method finds the right period only 
when all visible transits are overlayed, the answer is yes. Again, 
as in the case of a single observation, a requirement for the pe- 
riod detection is that the minimum of the global light curve is 
associated with a transit. 

We simulated three 150 d observations separated by gaps of 
the same length. The global light curve then lasts for about 2.5 
yr (see Figure |9]bef ore and after filtering). To ensure that at least 
one transit is visible in each data set, we set the first one at the 
half of the first observation and allowed the period to change 
over a low value P =300d±37.5d. In this exercice we only fo- 
cused on recover periods close to 150 d. 



Fig. 8. Same as Figure|7]but for low stellar activity. 

Figure[TO]and TablelH show the same results as Figure|2]and 
Table |2l respectively, but for multiple observations. The main 
conclusion from Figure [10] is that regular longlasting gaps do 
not affect the recovery of the transit period when close to 150 d. 
Results are similar to the case of one single data set with low 
stellar activity. The threshold for the minimum detected planet 
radius is again around 2 R^. 

6. Summary 

In this paper, we have presented a new method for detecting tran- 
sits of a planet orbiting around its parent star, and we investigate 
the limit of the smallest detectable planet when assuming that 
the star has similar characteristics to the Sun. We simulated light 
curves as described in Section|2l and, before applying the transit 
detection algorithm, we prefiltered the data to removing low and 
high frequency noise. 

Assuming a single observation of 150 d, the period of at least 
3 transits from a planet with a radius down to ~2 Re can be re- 
covered. The detection method does slightly depend on the ap- 
parent magnitude, but not on the transit duration in the given 
range of parameter values. The recovery of the right period does 
not depend on the number of transits, which only affects the con- 
fidence level (more transits will increase the confidence level). 

We have also finally investigated the case of multiple obser- 
vations (3xl50d) of the same field, when regular gaps of 150 d 
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Fig.9. Same as Figure[3]but in the case of three 150d observa- 
tions separated by gaps of 150 d. The phase of the first transit is 
at the half of the first data set. 
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Table 4. Same as Table|2]but in the case of three 150 d observa- 
tions separated by gaps of 150 d. 

Planet radius (R©) 
Radius NB. YES (%) NO (%) RA. (%) M.D. (%) 



1-2 


23 


9 


91 


4 





2-3 


23 


65 


35 


4 


4 


3-4 


29 


83 


17 


7 





4-5 


25 


96 


4 












Magi 


litude (Mag) 






Mag. 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


12-13 


28 


79 


21 








13-14 


24 


79 


21 





4 


14-15 


23 


65 


35 


4 





15-16 


25 


36 


64 


12 





Transit duration (hours) 


Tdur 


NB. 


YES (%) 


NO (%) 


RA. (%) 


M.D. (%) 


4-12 


32 


69 


31 


3 





12-20 


31 


65 


35 








20-28 


18 


50 


50 


11 





28-36 


19 


74 


26 


5 


5 
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Fig. 10. Same as Figs|2]and[8]but in the case of three 150d ob- 
servations separated by gaps of 150 d. 



are present and when one single transit is visible in all data sets. 
In this case, assuming the period is close to 150 d, the limit for 
the smallest detected planet does not differ from the case when 
all transits are observed in a single data set. These multiple ob- 
servations may represent an application for the CoRoT or simi- 
lar space missions to recover planetary transits when the orbital 
period is much longer than the single observation duration. 
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